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ABSTRACT 

We present spatially resolved imaging and integral field spectroscopy data for 450 cool giant 
stars within 1 pc from Sgr A*. We use the prominent CO bandheads to derive effective tempera- 
tures of individual giants. Additionally we present the deepest spectroscopic observation of the 
Galactic Center so far, probing the number of B9/A0 main sequence stars (2.2 — 2.8 Mq) in two 
deep fields. From spectro-photometry we construct a Hertzsprung-Russell diagram of the red gi- 
ant population and fit the observed diagram with model populations to derive the star formation 
history of the nuclear cluster. We find that (1) the average nuclear star-formation rate dropped 
from an initial maximum ~ 10 Gyrs ago to a deep minimum 1-2 Gyrs ago and increased again 
during the last few hundred Myrs, and (2) that roughly 80% of the stellar mass formed more 
than 5 Gyrs ago; (3) mass estimates within R ~ 1 pc from Sgr A* favor a dominant star formation 
mode with a 'normal' Chabrier/Kroupa initial mass function for the majority of the past star 
formation in the Galactic Center. The bulk stellar mass seems to have formed under conditions 
significantly different from the young stellar disks, perhaps because at the time of the formation 
of the nuclear cluster the massive black hole and its sphere of influence was much smaller than 
today. 



Subject headings: 
stars: early-type 
spectroscopy 



Galaxy: center — 
- stars: late-type 



galaxy: nuclear stellar cluster 
— stars: initial mass function - 



- galaxy: star formation — 
infrared: stars — infrared: 



Introduction 



The Milky Way nuclear star cluster (NSC) is of special interest since it is the closest galactic nucleus. 
It offers the unique possibility to resolve the stellar population and to study the composition and dynamics 
close to a central black hole at an unrivaled level of detail. The ability to resolve individual stars together 
with continuous mo nitoring of the innermos t stars has proven the existence of a 4.3 x 1O 6 M supermassive 
black hole (SMBH) (jEisenhauer et al.ll2005l : iGhez et al.ll2008t iGillessen et allEooi l beyond any reasonable 
doubt. Surveys of galaxies h ave found scaling rela tions between the bulge mass and the mass of the central 
massive object of a galaxy ( Ferrarese et al. 20061) . This massive object can be either a NSC or a SMBH, 
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depending on the bulge mass. SMBHs are typically found in massive bulges, while NSCs are common in 
low mass bulges. This indicates a mutual evolution of the bulge, the SMBH and the NSC. The Milky Way 
with its Bulge mass of 10 10 Mq falls right i n the transition region fr om galaxies, which are NSC dominated 
to galaxies , which are SMBH dominated ( Graham fe Spitler 20091) . For a cr itical discussion we ref er to 
Seth et all (|2008l , bold ). With a radius of 5 pc and a mass of 2 - 3 x 10 7 M (jLaunhardt et alj|2002h . the 



Milky Way NSC is typical. This lucky coincidence makes the Milky Way's NSC an ideal test case for nuclear 
co-evolution of galaxies in general. The evolution of a galactic bulge and the central massive object (NSC 
or SMBH) as well as the physical reason for the scaling relation is however poorly understood. When does a 
NSC form? Does it evolve simultaneously with the bulge or from a preexisting bulge? How does the SMBH 
influence t he star formation mod e? In fact, the age of a NSCs is the largest uncertainty in the determination 
of its mass ( Ferrarese et alj|2006[) . Deriving the star- formation history from integrated quantities as obtained 
from distant galaxies is highly model dependent and is especially uncertain for old ages. These limitations 
can be overcome in the GC because the stellar populations can be resolved. 

Numerous papers have studied the composition of the Milky Way's NSC. These studies have found that the 
stellar population can be divided into two classes. The cool and evolved giant stars and the hot and young 
main sequence /post-main sequence stars. The exis tence of massive young stars is evidence for very recent 
star formation ( Forrest et al Il987t Alle n et al . 1990). The most massive stars (WR/O stars) reside in a fairly 
complex structure, which may be described as a combination of a prominent warped d isk, a second disk-like 
structure high l y inclined relative to the main disk, and a more isotropic component ( Paumard et al. 20061 : 
Lu et alJliooa bartko et al.ll2009l ) at a projected distance of 0.8-12" from the SMBH Sgr A* (1" = 0.04 pc, 
assu ming a GC distance of 8.3 kpc). The GC disk features must have formed in a rapid star burst ~ 6 Myrs 
ago ( Paumard et al.lfeood : iBartko et al.ll2010 ). with a highly unusual initial mass function (IMF) that favored 
the formation of massive stars (dN/dm ~ m a ; a — —0.45). This extre me IMF deviates significantly from 
the standard Chabrier/Kroupa IMF with a powerlaw slope of a = —2.3 (|Kroupall200ll ; IChabrierl l2003') and 
seems to exist only in the vicinity of the S MBH. A less massive population of ordinary B-stars can be foun d 
in the innermost 1", the so called S-stars ( Eisenhauer et al. 2005 : Ghez et al. 20081 Gillessen et al. 20091 ) . 
The origin of the S-stars is a mystery, because the in-situ formation in the vicinity of the super massive black 
hole (SMBH) is very improbable. Yet, the K-band luminosity function (KLF) of the S-stars is consistent 
with a canonical Chabrier/Kroupa IMF. The largest population of the resolved stars however, are giants 
with masses between 1 and 2 M Q . Previous kinematic studies of this population have shown th at the cluster 



dynamics are consistent with a relaxed system, slowly rotating in the plane of the gal axy (Trippe et al 



20081 ISchodel et"aL 20091). Detailed ab undance determinations of luminous cool giants ( Carr et al. 2000; 
Ramirez et al.ll2000l : ICunha et al.l I2007T) found a metallicity distribution peaking at [Fe/H] = 0.14 ± 0.16, 



close to the solar value. 

The giant population is well suited for observations in the K-band. The giants have typical temperatures 
between 5100 and 2800 K. In this temperature range, the CO bandheads at 2.29 — 2.38 /im are the most 
prominent spectral lines i n the K-band. The equivalent widths of these lines correlates with temperature 
(jKleinmann fc Hall|[l986l I Wallace fc Hinklej|l996l iForster Schreiberlboool) . This allows the determination 
of individual temperatures for the red giants (RG). Together with photometric data for the luminosities 
it is thus possible to construct a Hertzsprung-Russell (H-R) diagram of the GC population and thus to 
constrain the star formation history in the immediate vicinity of the SMBH. Detailed s tudie s of the giant 



popul ation of the Galactic Center have been performed previously by iBlum et al.1 (|2003l ) and iManess et al 
(|2007l ). Using a magn itude limited samp le of 79 AGB and supergiant stars (50% complete at tuk < 10) 
within the central 5 pc. IBlum et al.l ([20031 ) found a variable SFR. They claimed that roughly 75% of the stars 
are older than 5Gyrs. An intermedi ate period of l o w sta r formation was followed by a recent (< 100 Myrs) 



period of increased star formation. IManess et al.1 ([20071 ) used 329 giant stars (ttik < 15.5), more than 5 
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magnitudes deeper than iBlum et al.l (|2003l ). albeit at the cost of covering only 12% of the total area in the 
central 1 pc. They found that the giant population is relatively warm, i.e young. Consequently they favored 
models with either a top-heavy IMF and a constant star formation rate, or with a canonical IMF and an 
inc reasing star formation r ate during the last few Gyrs. The top-heavy IMF scenario was recently challenged 
by Lockmann et al. ( 201ol) who argued that this would predict an overabundance of stellar remnants result- 
ing in a mass-to-light ratio larger than observed. 

The predicted density profile of a relaxed stellar system in the vicinity of a SMBH is a cusp with a radial 



density profile of n 



7 w 7/4 (Bahcall-Wolf cusp). Strictly speaking, this is only true for a single 



stellar mass cusp. In multiple mass configurations, 7 can range from 1/2 to 11/4. Recent observations have 
shown that the projected radial distribution of the giant p opulation is actually flat at small radii from Sgr A* 

There ma y even be a central hole in the 3D dis- 



|Po et al.ll2009HBuchholz et alj|2009t 



Bartko et al 



lant pc 
2010l). 



tribution of late type stars ( Buchholz et al. 20091 : Schodel et al. 20091) . In either case, the visible distribution 



seems quite inconsistent with a Bahcall-Wolf cusp. It is possible that a hidden Bahcall-Wolf cusp is present 
consisting of stellar remnants or unresolved main-sequence stars. However, this would require the stellar 
luminosity function to change in the few innermost arcseconds. The bright stars need to b e removed by 



mechanisms acting pr eferentially close to the SMBH. Collisional stripping of stellar envelopes (jGenzel et al 



1996; 



Alexander! 1999h is such a mechanism that keeps stars from r eaching their peak luminosity. Tidal strip- 
ping close to the SMBH can deplete the giants in a similar way ( Dayies et al.l 12006). Another explanation 
might be a top-heavy or truncated IMF ( Navakshin fc Sunvaev 20051) that produces preferentially massive 
stars as observed in the young disks. Massive s t ars q uickly become stellar remnants compared to evolved 
low mass stars that dominate the KLF. iMerrittl ([20091 ) argues that the absence of a Bahcall-Wolf cusp can 
be naturally explained if the stellar population is not relaxed. They find that the relaxation time scale in 
the innermost parsec is greater than 5 — 10 Gyrs. The most important question to discriminate between 
depletion processes and a relaxation time effect is therefore, whether the system is old enough to be relaxed. 



The paper is organized as follows; In section [2] we present the data and in section [3] the spectral 
classification and calibration of spectral indices. In section |4] we construct the H-R diagram that is then fit 
with model populations in section [SJ The star formation history, IMF and mass composition of the nuclear 
cluster is presented in section [6] The results are discussed and compared with other works in section [7] We 
conclude in section [8] 



2. Observations and data processing 

This work relies on spectroscopic and imaging data obtained at the VLT in Cerro Paranal Chile between 
2003 and 2010. The observations were carried out under the program-ids 075.B-0547, 076.B-0259, 077.B- 
0503, 179.B-0261 and 183.B-0100. 



2.1. Imaging and photometry 



The photometr ic data were obtained with the adaptive optics camera NACO ( Rousset et al. 20031 : 



Hartung et al.l 120031) . The photometric reference images were taken on the 29th of April 2006 and on the 



31st of March 2010. We used the H- and K'-band filter together wit h a pixel scale of 2 7 mas /pixel. To each 



image we applied sky-subtraction, bad-pixel and flat- field correction (jTrippe et al.ll2008l ) . All images of good 
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quality obtained during the same night were then combined to a mosaic with a FoV of 40 x 40". 



2.1.1. Diffuse background 



In terms of mass, the resolved stellar population represents only the tip of the iceberg. The bulk of the 
stellar mass is unresolved. Therefore the diffuse background emission of the GC contains valuable information 
on the cluster composition and its formation. Naturally it is very challenging to estimate the unresolved 
background in a crowded stellar field. A significant fraction of the background light originates from the 
uncorrected seeing halos of bright stars. Furthermore, anisoplanatism complicates precise photometry across 
the Fie ld of View (FoV ) when using AO. By using multiple point spread function (PSF) photometry on NACO 
images ISchodell d2010h tried to ov ercome these limitations. The resulting photometry for resolved stars was 
published in |§ch5del et all (fcoioh and a map containing the unresolved background flux was published in 
Schodel ( 2010l ). We used H-band data from both publications to derive the fraction of light contained in 
the diffuse background. The H-band has the advantage that it is least affected by the surrounding nebular 
emission. 



Since no completeness map was published bv lSchodell (120101) . we used a superb H-band mosaic with Strehl > 
20% from the 31st of March 2010. We applied StarF inder (Diolaiti et al. 2000l) and converted our detections 
in magnitudes by referencing them to Table A2 of ISchodel et al.l (|2010l ). To derive the completeness, we 
applied a common technique of creating and re-detecting artificial stars of various brightness. As a limiting 
magnitude, separating diffuse background from resolved sources, we used mH, cu t = 19.45 ± 0.12. The 
magnitude was chosen such that the total uncorrected flux of all detected stars is equal to the co mpleteness 
corre cted flux of stars brighter than the magnitude cut. Using an extinction of Ah = 4.65 ±0.12 (|Fritz et al 



2006 



|| Gillessen et al. 2009 ) we obtained an 



20111 ) and a GC distance of i? = 8.3 ± 0.35 kpc (|Ghez et al. 
absolute separating magnitude of Mjy,cut = 0.27 ± 0.34. The error estimate accounts for the absolute 
photometric uncertainty, the cut-off, extinction and distance uncertainties. We then subtracted the difference 
between the stellar flux (up to the limiting magnitude) in our H-band image and the stellar flux in the 
Schodel et all (l2010l) list from the diffuse background flux. Furthermore, we subtracted the known early- 
type stars since they dominate the light in the inner 10". 

We find that the diffuse background (m H > 19.45 ±0.12) contains H dm /H = 27 ±9% of the total (diffuse + 
resolved) H-band flux. The 1 a error contains the variation across the field (6%) and the uncertainty due to 
the cutoff between diffuse and resolved light Mjy jCUt (2%). Taking into account the completeness, we added 
another 3% systematic error due to undetected early-type stars. 



2.1.2. Mass-to-light ratio 



To derive the total H-band lumino s ity of the inner 1.2 pc, we used the resolved and unresolved stellar H- 
band flux stated in lSchodel et al.l (|2010l ); ISchodell (|2010l ). To account for the spatially varying extinction of the 
GC, we used their published extinction map. Unfortunately the extinction ma p is only available for K-band. 
Therefore we scaled the map to an average H-band extinction of Ah = 4.65 (jFritz et al.l 1201 II ). Correcting 
for extinction and taking the distance into account, we obtained a luminosity of H — 4.51 ± 0.54 x 10 6 Lh.q 
enclosed within a projected radius of r 2 D < 1.2 pc. From this we subtracted the contribution of known early- 
type stars (0.64 x 10 6 Lh,q). While the young stars are known to be close in 3D to the center, only 50±5% of 
the old stars at a projected distance of r 2 D < 1-2 pc are also contained within a 3D distance of 1.2 pc (assuming 
the radial density profile of Schodel et al.2007). Applying both corrections we find a total luminosity of 
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H(R < 1.2 pc) = 1.94 ± 0.33 x 10 6 L Q<H . Of this value 0.52 ± 0.2 x 10 6 L QM can be attributed to the diffuse 
background. The d ynamical mass enclosed (excluding the SMBH) is M(R < 1.2 pc) = 1.4 ± 0.7 x 10 6 M 
(jGenzel et al. 2010 '). Thus we derive a total mass-to-light ratio of M/H = 0.7± OAMq/Lq^h and a diffuse 
mass-to-light ratio of M/Ham = 2.6 ± UMq/Lq^. 

By assuming an average intrinsic colour of ttih-k = 0.2 we can convert the total H-band luminosity into 
a K-band luminosity of K(R < 1.2 pc) = 3.0 ± 0.44 x 10 6 L Q j<-. The corresponding mass-to-light ratio is 
M /K = 0.5±0.3 M /L QyK . The mass-to-diffuse-light ratio is M/K dm = 1.9±1.2 M Q /L Q ^ K . The bolometric 
mass-to-light ratio is M/Lboi = 0.7 ± 0.4 Mq/ Lq^ q \. 



2.2. Spectroscopy 



O ur spectroscopic data were obtained with th e adaptive optics assisted integral field spectrograph SIN- 
FONI (jEisenhauer et alj|200a bonnet et al.ll2004h . We used all high quality SINFONI data sets available to 
us, and taken between 2003 and 2010. In total we used 32 fields with pixel scales between 25 and 100 mas. 
The data output of SINFONI consists of cubes with two spatial axes and one spectral axis. Depending on 
the plate scale, an individual cube covers 0.8" x 0.8" or 3.2" x 3.2"; the spectral resolution varies between 
2000 and 4000 depending on the chosen bandpass and the field of view. Roughly 70% of t he stars used in 
this work were observed at a resolution of R ~ 2000. We used the data reduction SPRED ( Schreiber et al 



2004[|Abuter et allbood ). including bad-pixel correction, flat-fielding and sky subtraction. The wavelength 
scale was calibrated with emission line gas lamps and fine-tuned on the atmospheric OH lines. Finally we 
removed the atmospheric absorption features by dividing the spectra through a telluric spectrum obtained 
in the respective night. 



2.2.1. Deep spectroscopy of the GC 

Very good seeing conditions and long integration times on two SINFONI fields allowed us to perform 
the deepest census of the GC population so far. The two fields cover roughly 18 square arcseconds and are 
~ 50% complete down to rriK < 18. The location of the fields can be seen in Fig.[TJ Both fields were placed 
in regions with low extinction and without bright stars in the vicinity. One of the fields (East), at a distance 
of 7.4", probes the disk region. The second field (North) is located at the outer rim of the disks at 13.5". The 
richness of faint stars can be seen in Fig.[TJ The identification of a rax < 17 star as late-type requires only 
two hours of integration on source. The prominent CO bandheads are easy to identify and the multiplicity 
of the lines helps to rule out misidentifications in noisy spectra. However, the identification of an early-type 
star is much more challenging. The early-type stars are identified by the hydrogen Br7 line. Unfortunately 
the line is less prominent than the CO bandheads. The identification is further hindered by the surrounding 
nebular Br7 emission. This emission is patchy and varies significantly in intensity and velocity on scales 
of 1". Depending on the location of the star, this can lead to a wrong background subtraction during the 
data analysis mocking a stellar line. Because the stellar absorption line in A-stars is a broad Lorentzian, 
it is clearly different from the narrow Gaussian emission line of the gaseous background. Thus in case of 
sufficient SNR an identification is possible. To achieve a sufficient signal-to-noise ratio (SNR) an integration 
time of about five hours on source is required to identify an early-type star as such with mx < 18. 
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Fig. 1. — Three-color image of the innermost 1 pc of the Milky Way (K-band: blue; L-band: red). The white 
cross indicates the position of Sgr A*. The U-shaped nebular emission visible in the L-band, the so called 
Mini-Spiral, is powered by the ionizing radiation of the young O/WR stars in the vicinity. The position of the 
two deep spectroscopic fields is indicated. All photometrically detected stars not confused with neighboring 
stars are indicated by open circles. The color coding resembles the stellar type. Cold giants with confirmed 
CO bandheads are indicated in red. Hot main sequence stars identified by Br7 absorption are indicated in 
blue. Green stars are not identified. 
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2.3. Source selection 



The selection of stars used in the H-R diagram fitting is based o n a master list of ~ 6000 stars at 
projected radii between 0.1" < r < 25" (or 4 x 10~ 3 to 1 pc), for which Trippe et al. ( 20081) derived proper 



motions as well as K-band photometry. The list consists of well isolated stars identified on several images 
without overlapping neighbors (separation > 130 mas). For our final analysis we retained only those entries 
in the master list that also had an unambiguous identification in a corresponding SINFONI cube. In total 
we collected ~ 1300 spectra. After removing duplicates (due to overlapping fields), roughly 1000 spectra 
remained. Out of this sample we used only stars with a SNR > 10 and a CO index EW(CO) > 3.5 (i.e. 
no early- type stars). The SNR was measured in the continuum bands stated in Table[TJ This left us with 
about 800 giants for which we were able to determine individual temperatures and velocities. In order to 
have a homogeneous and deep sample, we used in the H-R diagram fit only stars contained in fields with 
a completeness > 50% at vtlk = 15 (see next section). The remaining sample consisted of 450 stars with 
magnitudes between 0.5 < Mk < — 8. The contribution of fore- and background stars is negligible in the 



GC. The contamination due to stars which are not members of the cluster is of the order 1% (jBuchholz et al 
2009h . 



2.4. Detection probability 

The spectroscopic fields used in this work differ in spatial resolution, covered area and integration time. 
The main limitation for spectroscopy in the GC is stellar crowding. This is why we probed most of the 
inner lpc with the small scale (100, 25 mas) AO-assisted modes of SINFONI. For these observations, we 
used typically one hour of integration time on source. Two additional cubes (100 mas) are exceptionally 
deep with integration times of 4.5 and 7.8 hours on source. The spectroscopic completeness of these cubes 
is > 50% at itlk < 18. The results of the deep observations are discussed in Sec. 13.11 The spectroscopic 
completeness was determined by comparing the total number of stars contained in the master list, with 
the number of stars for which we could extract spectra. The photometric completeness of the images was 
determined as described in Sec. 12.1.1] All observed fields and the combined completeness (for itlk = 15) 
are shown in Fig. [2] Fields with a high completeness of up to 80% are found at several arcsecond distance 
from SgrA*. Further in, the completeness degrades due to stellar crowding. Only the inner one arcsecond 
is sampled with the smallest pixel scale (25 mas) and reaches a completeness of about 50%. 



Table 1. Equivalent width measurement intervals 



Feature 


Wavelengths (/im) 


12 CO(2,0) band 


2.2910-2.3020 


12 CO(2,0) continuum I 


2.2300-2.2370 


12 CO(2,0) continuum II 


2.2420-2.2580 


12 CO(2,0) continuum III 


2.2680-2.2790 


12 CO(2,0) continuum IV 


2.2840-2.2910 
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RA foresee] 



Fig. 2. — The AO assisted SINFONI fields overlaid on a GC image. The color of the fields indicates the 
combined spectroscopic and photometric completeness for a ttlk = 15 star. The completeness is mainly 
driven by the exposure time of the individual fields and the observing conditions of the respective night. 
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3. Spectral classification 

The K-band provides two prominent spectral features allowing a spectral classification. Early-type stars 
(T > 5000 K) can be identified by the presence of a Br7 absorption line. Late-type stars (T < 5000 K) on 
the other hand show weak or no Br7 absorption but strong absorption line blends known as CO bandheads. 
The bandheads appear at the spectral type ~ G4 and increase in strength up to M7. The strength also 
increases with luminosity class from dwarfs to giants and to supergiants. Yet, for this work only giants are 
of interest. Main sequence stars with spectral type G4 and cooler are fainter (ttik > 20) than the detection 
limit. S upergiants of the s ame spectral type on the other hand are extremely bright (ma < 9) and easy to 
identify |Blum et al.ll2003h . 



3.1. Deep census of the GC population 

The two deep S1NFONI fields allowed a deep census of the GC population with 50% completeness down 
to ttik < 18. 



3.1.1. A-star detection 



Fig. Q] shows the deep fields. Red and blue circles mark stars identified as late-type or early-type. One 
star was not confirmed as early- type (green circle), but showed no CO bandheads and is therefore deemed to 
be an early- type candidate. Within the 18 square arcseconds of the two fields, we thus detected five early- 
type stars and on e cand idate. Their spectral type was determined by comparis on with temp l ate sp ectra of 
Wallace fe Hinkld (1997) and the absolute K-band magnitude Mjf taken from lAllen fc Coxl ( 20001 ). Three 



of them were brighter than < 17 and were identified as dwarfs with spectral types between B2 and B8. 
The two faintest early type stars identified were B9/A0 dwarfs with rriK — 17.2 and — 17.6. These 
stars are the faintest main-sequence stars reported so far in the Galactic Center. Stars of that spectral 
type have main-sequence lifetimes between 360-730 Myrs and masses between 2.2 — 2.8 M© (jAllen fe Cox 
20001 ). The spectra of the two A-stars are shown in Fig.[3] For comparison, two spectra of late- type giants of 
equal K-band luminosity are shown. The candidate is of similar luminosity falls in the same main-sequence 
category. 



3.1.2. K-band luminosity function 



The deep census of the GC allowed us to construct a KLF with 12 < tjik < 18 (Fig. [4}. The best fitting 
slope of the lumin osity distribution d log N/dm k — (3 is found to be (3 — 0.33 ± 0.03. This is consistent 
with the results of iBuchholz et al.l (120091) . who determined th e slope of stars with a limit i ng magnitude of 
nig < 15.5. As already discussed by Alexander fc Sternberg ( 1999 ); Genzel et al. ( 20031) : Buchholz et al. 
(|2009l ). the slope is typical for an old, bulge- like population. This is in good agreement with our H-R 



diagram fitting (Sec. l5.lj ) showing that most of the star formation occurred more than 5 Gyrs ago. However, 
the KLF alone is not sufficient to constrain the star formation history for an old population. The KLF 
alone cannot cons train the actual forma tion history and the IMF since the KLF slope is very insensitive to 
both parameters ( Lockmann et al. 2010l) . For example, the admixture of young and bright giants is hardly 
detectable in the KLF due to their low number. Yet, they can be easily identified in the H-R diagram. 
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Fig. 3. — Spectra of the faint A0/B9V main-sequence stars in th e Galactic Center (uppe r panels). The red 
line shows B9/A0 V template spectra of the solar neighborhood ( Wallace fc Hinklel 1997 ). The dashed bars 
indicate the width of the nebular emission in the vicinity of the stars. The name in the upper left corner 
relates the spectrum to the star in Fig. [I] For comparison, two late-type spectra of similar brightness ares 
shown in the lower panels. The corresponding late-type templates are indicated in red. 
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Table 2: Early type stars contained in the deep spectroscopic fields. 



Field 


mK 


RA a ["] 


Dcc a ["] 


v x [mas/yr] 


v y [mas/yr] 


v z [km/s] 


Disk member b 


North 


14.1 


0.0 


12.3 


2.44 ±0.10 


0.55 ±0.24 


-100 ±40 


CCWS 


North 


17.2 


2.0 


13.8 


0.62 ±0.27 


0.96 ±0.4 


116 ±55 


no 


East 


14.2 


6.5 


-2.7 


-2.73 ±0.14 


-3.73 ±0.26 


-60 ± 7 


no 


East 


16.4 


7.3 


-4.1 


2.89 ±0.23 


-1.25 ±0.16 


-163 ±44 


no 


East 


17.7 


8.0 


-2.7 


2.19 ±0.39 


-0.12 ±0.24 


-41 ± 52 


no 



"Distance in arcseconds relative to Sgr A*. 

^Potential members of the clockwise disk (CWS) or the counter-clockwise disk (CCWS) as defined bv lBartko et all l|2009l ). 



3.1.3. Origin: Disk members? 

One of the observed fields is centered on the disk region, while the second field probes the outer rim 
of the disks. Thus the early-type stars in the two observ ed fields are potent ial disk members. To test 
for membership, we used the same method as described in iBartko et al.l ( 201oh . The 3D velocity and the 
positions of all five stars are given in TableO The last column indicates a likely disk membership. We find 
that the two A-stars are clearly not disk members. Among the b righter B-stars , only one is consistent with 
the counter-clockwise disk (CCWS). This supports the results of IBartko et al.l (|2010h who find a disk IMF 
skewed towards massive stars. Their disk KLF predicts roughly the same number of stars with ttik = 14 
as with m,K — 17. A Chabricr/Kroupa IMF on the other hand predicts five times more mx — 17 stars. 
Although the statistical significance of only five early-type stars is limited, it is intriguing that we find one 
rtiK — 14 disk member (consistent with the average B-star de nsity) but no fainter disk members. The faint 
B/ A-stars exhibit an isotropic orientation ( Bartko et al. 2010f) and are most likely the remains of older star 
bursts. 



3.2. Supergiant IRS7 

Two red supergiants are known to reside within 1 pc of Sgr A* (jBlum et al.ll2003l: baumard et al.lbooll) . 
Those supergiants are significantly younger than the rest of the old ste llar population. The supergiant IRS 7 

Carr et al.1 (|2000h find an initial mass of 20 M Q 



is the youngest of this particular luminosity class. iCarr et al.1 (|2000l ) find an initial mass of 20 M Q and 
stellar age of ~ 7Myrs and attribute it to the recent starburst. This is supported by its kinematics (Bartko; 
pri v. comm.) placing I RS 7 on the clockwise d isk system. We simulated a 6 Myr star burst as found 



bv IPaumard et al.l |2006l) and IBartko et al.l |2010l) with the IAC-STAR code (see following sections). The 
predicted ratio of red supergiants (T c ff < 5000 K; ttik < 9) to blue supergiants (T ff > 6000 K; 9 < ttik < 14) 
is between 0.4-0.5% (independent of the assumed IMF). This is in good agreement with the observed ratio of 
118 blue supergiants with rtiK < 14 ( Bartko et al. 201oTl to one supergiant, i.e. a ratio of 0.9%. The second 



supergiant has an age of a few ten Myrs ( Blum et al.l 



20031 ) and is thus older than the disks. 
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3.3. CO index definition 



The 12 CO(2,0) bandhead has been widely used as a temperature indicator. This has brought up nu- 
merous defini tions of CO indices. A detaile d comparison of regularly used index definitions was recently 
performed by Marmol-Queralto et al. ( 20081) . The analysis showed that some index definitions are system- 
atically more affected by spectral resolution, velocity error, curvature of the spectrum and SNR than other 
definitions. In particular, the index definition of Klcinman & Hall 1986 proved to be very sensitive to those 
effects. The index measures the continuum and the line flux in two narr ow bands in close p roximit y.A similar 
i ndex was used in previous studies of the Galactic Center po pulation by Maness et al. (|2007l ) and I Blum et al 



(2003). However the authors used a wider bandpass than iKleinmann fc Halll (1986) (0.015 /j,m instead of 
0.0052/im). In the following the index is referred to as (BL03). Systematic errors in the measurement of the 
CO strength lead to systematic errors in the temperature estimation. Since the temperature is a tracer of the 
stellar age, this might lead to a bias in the determinatio n of the star format ion history. Indeed, the paucity 
of cool stars in the GC is the chief constraint which led iManess et al.l ([20071 ) to conclude a top-heavy IMF 
may have persisted throughout the star formation history of the Galactic Center. To estimate the impact 
of systematics in light of this recent work, we tested the index BL03 against various effects and compared 
its performance with the alternate index proposed by Frogel et al. (2001, FR01). It uses several narrow 
bandpasses to estimate the CO continuum with a linear fit (see Tab. [I}. The line and continuum regions for 
both index definitions can be seen in Fig. [5] 



3.3.1. Index computation 



Before computing the indices, each stellar spectrum was shifted to rest-wavelength. We then divided 
the spectrum by a second-order polynomial fit to remove the curvature of the spectrum. For the continuum 
fit, we excluded regions with significant absorption lines; the polynomial fit is necessary to account for the 
the large (Ak ~ 3) and spatially v ariable extinction in the GC. We then computed the BL03 index according 
to the recipe of 



Blum et al 



(2003J). The FR01 index was computed in a similar way, with the bandpasses as 



described in iFrogel et al.l (|2001[ ). The continuum level w c (X) is estimated with a linear fit to the intervals 
stated in Table[TJ The equivalent width is measured according to: 



EW(CO) 



band 



WC ~ Wn nc 
U> C 



)dA. 



(1) 



We tested both indices (BL03, FR01) with template spectra of different resolution. We also reddened the 
template spectra artificially to determine the impact of extinction. 



3.3.2. Systematic error sources 

Reducing the resolution from R ~ 3000 to R ~ 2000 decreases the BL03 CO index by a factor of « 0.92. 
The artificial reddening and the corresponding change of curvature of the spectrum, causes a reduction of 
« 0.94. Both effects cause the BL03 CO index to be underestimated by a factor between 0.90-0.85. The 
systematic underestimation of the BL03 CO strength leads to an overestimation of the stellar temperatures 
by ~ 200 K. The FR01 index decreased by less than a factor 0.98 due to extinction. Degrading the resolution 
from R ~ 3000 to R ~ 2000 shows no measurable impact. We estimate the combined systematic effect to 
be less than 0.97 (50 K). The BL03 index suffers also from contamination of the Mgl line contained in 
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Fig. 4. — Completeness corrected KLF of the confirmed early- (blue) and late-type (red) stars found within 
the two deep spectroscopic fields. The green KLF shows the early type stars including the unidentified star. 
The solid line indicates the best-fit powerlaw to the late- type KLF. 
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Fig. 5. — Left: Index definition of lFrogel et al.l (|200ll ). The gray area indicates th e line region and the whit e 
area indicates the continuum regions use d for the fit . Righ t: In dex definition of Kleinmann fc Halll (|1986). 
An adaptation of that index was used in iBlum et al. | |2003l) and lManess et al.l |2007t) . 
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the continuum bandpass of the index. The average Mgl line strength of the GC giants is similar to the 
calibration giants. Thus, the continuum estimation is not biased. Yet, the line shows some intrinsic scatter 
introducing a statistical error in the CO index computation. To minimize the statistical and syst ematic error 



source s, we therefore adopted the index definition of FR01. We have to note that the results of I Blum et al 



(|2003l ) did not suffer from the resolution dependence of their CO index because the bulk of their GC spectra 
were of the same resolution as the comparison stars. 



3.4. Temperature calibration 



For the temperature calibration we used stellar spectra of giants with known T Q g. In total we used 33 
giants with spectra l types G0-M7 and meta llicities —0.3 < [Fe/H] < 0.2. The spectra were obt ained from the 
NOA O IR-libra rv dWallace fc Hinklej|l997h with R - 3000; the IRTF-library with R - 2000 (jRavner et al. 
2009) and from iForster Schreiberl |2000j) with R ~ 2000. The corresponding temperatures were obtained 



with the help of the SIMBAD database. The CO-T e fj relation for the template giants can be seen in Fig. [5] 
The best-fit to the data in the range 3.5 < CO < 24 is (where CO is the FR01 index in units of [A]): 



T cS = 5832 ±13 ° - 208. 25 ±32 84 • CO + 11.38 ±248 ■ CO 2 - 0.34 



±0.06 



CO J 



(2) 



The residual scatter is 119 K. The average metallicity of the calibration star s of [Fe/H] = —0.1 is somewhat 
below the average metallicity in the GC of [Fc/H] = 0.14 (jCunha et al.ll2007l ). To account for the metallicity 
dependence of the CO stre ngth, we add a constant offse t of 50 K to the temperatures of the GC giants 
(estimated with the help of Marmol-Queralto et al. ( 20081 ); their Fig. 14). 



3.4-1- Systematic temperature uncertainty 

The systematic errors in the derivation of the temperature are the main drivers in the age uncertainty of 
a giant population (at a given metallicity). As discussed in the previous section, the CO index computation 
might be biased on a ~ 50 K level. More important however is the theoretical uncertainty. The Padua and 
Geneva Isochrones can differ by up to 80 K for the same star. Strictly speaking, this is not an observational 
uncertainty but in the fitting procedure the theoretical uncertainty must be taken into account. We included 
the bias by assuming a total systematic uncertainty of 100 K. To visualize the impact of temperature uncer- 
tainties, Fig.[7]shows the theoretical age vs. median temperature relation for red clump stars (0 < Mboi < 1)- 
Especially for old ages, the age interpretation is very sensitive to temperature changes. Small temperature 
biases of the order 100 K might change the derived age by several Gyrs. In the star formation history fitting, 
we used not only red clump stars. However red clump stars make up more than one third of all observed 
giants. Thus, the systematic uncertainty of the red clump population is representative for the whole giant 
population. 



3.5. Red clump spectrum: evidence for old ages and near solar metallicity 

The majority of the GC stars are red giants (RG). The RG phase is associated with the hydrogen-shell 
burning phase of a star. While the burning shell propagates outward through the stellar interior, the star 
expands and cools down. During this evolution, the star increases its brightness and climbs up the RG 
branch. Discontinuities in the stellar interior cause the brightness of the star to drop temporarily, when the 
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Fig. 6. — The EW(CO)-T e g- relation for the template giants is shown. Red triangles indicate spectra with a 
resolution of R ~ 3000, black dots represent spectra with a resolution of R ~ 2000. One 10 a outlier can be 
seen. This particular star has a metallicity ([Fe/H] = —0.5) significantly lower than the other template stars 
(< [Fe/H] >= —0.1). It was only included for comparison, but was not considered in the fit. The solid line 
shows the best fit to the data. The fit is reliable between 3 < EW(CO) < 24. The fit residual RMS is 119 K. 
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Fig. 7. — Relation between age and median temperature of the Red Clump is shown (solid line). Overplotted 
is the relation between age and mean initial mass of the Red Clump population (dashed line). The shaded 
area indicates the T e e uncertainty. The plot illustrates the temperature sensitivity to the age of the red 
clump population. For ages >5 Gyrs, the age resolution decreases significantly. 



burning shell passes a discontinuity. This happens several times during the ascent of the RG branch. This 
creates an overdensity of stars at a typical luminosity referred to as the red clump (RC). 



3.5.1. GC red clump 



The observed K-band magnitude and the spectral CO index allowed the determination of the bolometric 
magnitude and temperature for individual GC stars. The bolometric magnitude Afboi follows directly from 
the intrinsic, de-reddened Mr-. To derive the intrinsic Mk w e used a distance mo dulus of 14.6 (8.3 kpc). The 
extinction was calculated for each star individually, following Bartko et al. ( 2010l) . We used for each star the 
20 nearest neighbors to derive the median H-K colour. By assuming an intrinsic colour of Ams-K = 0.1 
(typical for K giants) we derived the reddening. By assuming a powerlaw extinction A\ ~ A~ 2,21 (Schodel et 
al. 2010, Fritz et al.(2011) we com puted the Ak extinction for each star. The bolometric correction (BCk) 



Blum et all |2003l) \BC k = 2.6 - (T Kff -_3800)/1500]. The assumed BC K is consistent with 



was taken from 

the bolometric correction library from iLeieune et al] ([19971 ) used in the subsequent population synthesis 
models. The bolometric 1 a uncertainty is 0.32 magnitudes. It is caused by the photometric error, the 
temperature- and the extinction uncertainty. The stellar T e g follows directly from the CO strength. The 
statistical 1 a T e g uncertainty derived by error propagation of the measured CO index uncertainty is 200 K. 
The systematic T e s uncertainty is 100 K (see Sec. 13. 4. J|) . The red clump of the GC giants can be found at a 
luminosity of Mbol = 0.6 and a temperature of T e g = 4310 K (see Fig. fTU)) . Since the red clump is narrow, 
both in magnitude and temperature, i.e. all stars have very similar spectral features, it is a viable approach 
to construct a median spectrum of the GC red clump stars. In order to construct a high SNR spectrum, we 
used 234 red clump stars with bolometric magnitudes < Mbol < 1 corresponding to 15.1 < uik < 16.1. 
The average SNR of individual spectra is between 10 and 15. All spectra were normalized and shifted to 
rest-wavelength as described in the previous section. Most of the spectra have a resolution of R ~ 2000. 
Spectra with a higher resolution were interpolated to match the latter. The equivalent integration time of 
the combined spectrum is of the order 300 hours. Figure|8] shows the resulting median spectrum of the GC 
giants. 



3.5.2. Temperature estimate from high SNR spectrum 

We compared the median spectrum with available library stars. The template star matching the co- 
added spectrum most closely is a K3III giant with a temperature of 4330 K. The agreement of the spectra 
is impressive. Only the sharpest peaks are slightly smeared out due to a residual velocity error. The 
12 CO bandheads at 2.294, 2.322, 2.352 and 2.383 /xm are well represented. Even the weaker bandheads of 
13 CO at 2.345 and 2.373 /im match the template (Fig. [5]). The comparison with colder and warmer giants 
(Fig. [5]) shows the sensitivity of the CO bandheads to the stellar temperature. The median temperature of 
the red clump as derived by the stacked spectrum T g = 4330 K agrees very well with the median value 
(T e g = 4310 K) of the individual low SNR spectra. The Cal and Nal line blends are moderately sensitive to 
temperature as well as to the surface gravity log<?. The Cal lines are in good agreement with the template. 
The same is true for the iron line blends at 2.22 and 2.23 /mi as well as Mgl, Si I and All. Unfortunately 
almost all of the atomic lines are contaminated with weaker atomic lines and lines of CN. This excludes the 
possibility of a spectral synthesis fit at the given resolution. However qualitatively, it is clear that the GC 



spectr um is closely matched by a solar metallicity spectrum, in agreement with the work of ICunha et al 



(120071) 
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Fig. 8. — The median spectrum of 234 GC red clump giants (black). Each individual spectrum has been 
normalized and shifted to rest-wavelength. The median spectrum of the GC red clump is closely matched 
by a K3 giant of the solar neighborhood (red) with a temperature of 4330 K. The most prominent lines are 
indicated. Almost every line is actually a blend of atomic and molecular lines. In case of a line blend, the 
most important line is indicated. 




Fig. 9. — The temperature dependence of the CO bandheads is shown. The median spectrum of 234 GC 
red clump giants (blue) is compared with three giants of the solar neighborhood with temperatures between 
4480 and 3900 K (gray shaded areas). 
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3. 5. 3. Sodium line blend: Potential evidence for mixing processes 



The strongest deviation between our average spectrum and the K3III template is found in the Nal 
lines. They seem to be intrinsically stronger than in the solar neighborhood. None of the template spectra 
in the temperature range of the red clump can rep roduce the Nal st rength. An increased Nal strength 
was reported previously in low- resolution spectra bv lBlum et al.1 (J1996). The Nal lines are actually blends 
of a couple of atomic lines and CN lines. The individual contribution depends on the stellar temperature. 
For temperatures similar to a K3III giant, sodium is the most important line. Cooler spectra are heavily 
influenced by Sc and the cool e st spe ct ra are domi n ated by CN. For a detailed analysis of lines in K-band 



spectra see Wallace fc Hinklel ( 19961) . Carr et al. ( 2000l) used high resolution spectra of IRS 7 to derive 



abundances. They found that the increased line strength of the Nal complex is mainly caused by stronger 
CN lines. The CN lines reflect extreme CNO abundances in the atmosphere of IRS 7, which they claim is 
probably the result of increased rotational mixing. 

Cunha et al.l (|2007T) also find CN-cycled material in the atmospheres of three luminous GC giants, however 
not as deeply mixed as in IRS 7. Increased rotational mixing, is predicted for de nse stellar cluster s, where tidal 
spin- up can lead to significantly higher rotation speeds of main-sequence stars (| Alexanderll2005l ) . This might 
explain the increased CNO material in the outer atmosphere layers of the luminous GC giants. Evidence for 
increased rotational mixing taking place in a supergiant like IRS 7 with a mass of 20 Mq and a lifetime of 
a few Myrs, are hard to transfer to the red clump giants of the GC with masses of only 1-2 M Q and ages of 
Gyrs. Th is is especially true since IRS 7 is significantly cooler (3600 K) than the red clump. 
However. I Maeder fc Mevnet (|2000h note that red giants with masses M < 1.5 Mq are susceptible to extra- 
mixing and lAlexanderl (|2005l) notes that tidal spin-up is most effective in long-lived low mass stars. Thus, 
the strong sodium lines might indicate that the red clump stars have undergone a mixing process different 
from the solar neighborhood, for example due to tidal spin-up. In principle, the sodium line strength can 
also reflect a peculiar chemical composition of the Galactic Center. With the current knowledge we cannot 
rule out that possibility. However, rotational mixing provides a convincing explanation for the increased 
sodium strength. This might provide interesting conclusions on the stellar evolution since fast rotating stars 
tend to be more luminous and redder and c an provide evide nce for the existence of an underlying dense 
stellar cusp of low-mass main-sequence stars (|Alexanderll2005l ). 

Apart from the Nal lines, the average spectrum shows no peculiarities. Small deviations of the GC spectrum 
at 2.219 and 2.349 /im are caused by poorly subtracted nebular emission of the Mini-Spiral. The deviation at 
2.317 /xm on the other hand coincides with a telluric line. Overall the agreement is intriguing and validates 
the method of constructing a median spectrum a-posteriori. 



4. Construction of the H-R diagram 



We constructed an H-R diagram for the GC cluster using the same distance modulus and extinction as 
described in Sec.13.5.11 The H-R diagram is displayed in Fig.[TUJ For comparison, we also included the data of 
Blum et all I^OO^ probing the luminous giants and supergiants in the i nner 2.5 pc. The co mparison data were 
rescaled to the updated extinction values used in this work. Although iBlum et al.l (|2003[) and this work rely 
on different te mperature calibrati ons, the overlap regions match very well. The GC population is compared 
to isochrones (|Bertelli et al.lll994l ) with Z = Z& and Z = 2.5 Z Q metallicity (Fig.[10l left). The right panel 
of Fig. HOI shows the comparison with with a synthetic model population. The model p opulation represent s 
the case of continuous star formation over 12 Gyrs with a GC metallicity of Z = 1.5Z & (jCunha et alj |2007). 
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Fig. 10. — Left: H-R diagram of the cool Galactic Cen ter stars; this work (red crosses). Stars younger than 
0.5 Gyr are indicated by star symbols; data points from lBlum et al.l (|2003l ) are indicated as black diamonds. 
Overplotted are isochrones for ages (from left to right) 0.05, 0.1, 0.5, 2 and 12Gyrs with solar metallicity 
(black, solid) and 2.5 times solar metallicity (blue, dashed). The upper black cross shows the statistical error 
and the lower cross shows the systematic uncertainty. Right: Observed H-R diagram (red crosses) overlaid 
on a synthetic mode l population assuming 12 Gyrs of continuous star formation; the metallicity is Z — 1.5 Zq 
(jCunha et al.ll2007D ; Different age intervals of the synthetic population are indicated; 5 < t < 12 Gyrs (dark 
gray); 1 < t < 5 Gyrs (green); 0.2 < t < 1 Gyrs (gr a y); 0. 05 < t < 0.2 Gyrs (blu e). The synthe t ic pop ulation 



was created with the co de of | A paricio fc Gallartl (|2004). The code used the iBertelli et al 
evolution library and the Leieune et al. ( 19971) bolometric correction library. 



(19941) stellar 
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4.1. Features of the H-R diagram 

The H-R diagram shows several distinct features. The most prominent feature is the RG branch con- 
sisting of old (> 1 Gyr) stars with masses of the order 1M . The GC red clump can be identified as an 
overdensity of stars at a lu minosity of -M h oi = 0.6. This agrees very well with the red clump in the solar 



neighborhood measured bv iGroenewegenl (|2008l ). The mean red clump temperature derived by individual 
spectra agrees very well with the estimate based on the high SNR median red clump spectrum (Fig. l3.5p . 
A second overdensity can be found at M\, \ = —0.8 and T c s — 3900 K. This feature is sometimes referred 
to as the AGB bump. All afore mentioned features are tracers of an old population. The H-R diagram 
shows a second branch of giants at T c g — 4800 K. This warm giant population is bright (Afboi < 0) and 
separated from the cold (old) branch in the same magnitude range. Yet, only ~ 10% of the giants between 
< Mboi < —4 (corresponding to 11.5 < tuk < 15.5) can be attributed to the young branch. Similar 
features are known in globular clusters with multiple stellar populations. The comparison of the GC H-R 
diagram with Chabrier/Kroupa populations of different ages can be seen in Fig. [TU] (right). Several obvious 
features are matched. The red clump is well represented in the data. The cool branch of giants with ages 
greater than 1 Gyr is obvious. The temperature and luminosity of the warm branch is matched by giants 
younger than 0.5 Gyrs. A deficiency of stars with ages of ~1 Gyr can be seen as a gap in the diagram. We 
don't detect a significant Horizontal Branch of old and metal poor giants as would be typical for old globular 
clusters. This supports the assumption of a predominantly old and metal rich population. 



4- 1.1- Outliers 



The coolest temperatures we find cannot be reproduced by even the oldest solar metallicity isochrones. 
Those outliers might be explained by a population of metal rich stars (Z > 2.5 Zq). However, this does not 
contradict the current assumption of a near solar metallicity in the GC. The cold outliers account for less 
than 8% of the total population. Since metallicity studies in the GC (jRarmrez et allbood Iciinha et al.ll2007h 
have probed not more than ten late-type stars, a small metal rich population might have gone undetected. 
This can be compared to the situation in the Bulge (Baade's Window), where a high metallicity tail of 
stars [Z > 2.5 Zq) makes up about 15% of the population (jZoccali et al.l 120081 ). Another explanation 
might involve stellar model uncertainties. Most of the outliers show the same luminosity as bright and cool 
AGB stars. Stellar models of these stars still suffer from large uncertainties. Model isochrones treat this 
evolutionary phase in a simplified manner and are thus rather unreliable. Furthermore, stars in these stages 
are known to pulsate with perio ds of hundreds of days. Du ring the pulsation, an individual star can change 
its temperature by up to 500 K (jLancon fc Mouhcinell2002l) . All systematic effects considered by us, lead to 
an overestimation rather than an underestimation of the stellar temperature. Yet, some of the stars can be 
heavily dust obscured and thus appear too faint for the given temperature. In the following fitting procedure 
we ignored stars that were cooler than allowed by the isochrones and the temperature uncertainty. This 
excludes only a small number of the old stars. Therefore we are confident that the impact on the results is 
negligible. 



4-1-2. The young giant branch 



Roughly 10% of the observed GC giants with magnitudes between < Afboi < —4 appear to belong 
to a branch of young (< 500Myrs) giants. With masses between 2.5 M Q < M < 6 Mq those stars are 
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descendants of main sequence B-stars. They are tracers of an intermediate age population in the Galactic 
Center. Their stellar age is smaller than the t ypical non-resonant two-body relaxation time in the Galactic 
Center of the order t m > few Gyrs in the GC |Alexanderll2005t lMerrittll2009l) . However, close to the SMBH, 
orbits are near-Keplerian. This causes interactions between stars to build up coherently. The randomiza- 
tion of the angular momentum vectors happens therefore on the fast vector resonant-relaxation timescale 



t V r < few ten Myrs (jHopman fc Alexander 



2006). 



However, the eccentricity and semi- majo r axis distribution is randomized on the slow scalar resonant 
timescale t SI (jHopman fc Alexander! 120061 ) . comparable to £ nr for r > 0.1 pc. The mean stellar mass of 
the young giants (3.2 M©) is significantly larger than the average mass of the old giants (1 M Q ). The more 
massive stars tend to sink to the center of a cluster due to dynamical friction exerted through the drag of 
lighter background objects. The mass segregation timescale in the Galactic Center is larger than in normal 
stellar clusters due to the presence of the SMBH and the corresponding higher velocity dispersion. The 
mass segregation timescale t s ~ t nT j4- scales with the relaxation time t nl but is also a function of the in- 
dividual stellar mass M* and the mean stellar mass M ([Alexander! 120051 1. Assuming the mean stellar mass 
is M ps IMq, the segregation time is t s ±i„ r > 1 Gyr. This is still significantly larger than the age of 
the young giants. Thus the radial distribution could not have changed significantly since the stars formed 
(neither through relaxation, nor through mass segregation). Yet, the angular momentum vectors (i.e the 
orientation of the orbits) will have undergone randomization. 



4-1-3. Kinematics and distribution 
To assess the dynamica l state of the young giants we computed their orbital distribution. We followed 



the method of iBartko et al.l (|2009i ). who used statistical arguments to infer the 3D distribution of a stellar 
population out of the 3D velocity and projected distance distribution. A coherent motion within a stellar 
system can be detected as a statistical preference for a angular momentum direction. However the orbital 
distribution of the young giants shows no significant excess as would be the case for a disk or a streamer. 
The orientation of the angular momentum vectors is consistent with being isotropic, as expected due to 
the fast vector resonant relaxation. The radial distribution of the young giants compared to the cool giant 
population is shown in Fig. 111! 

Both populations exhibit the same radial density distribution. The d ensity distribution of the total giant 
popu l ation has been descr ibed as a core or even a hole close to the SMBH (jBuchholz et al.ll2009t IBartko et al 



2010t ISchodel et al.1 120101 ) . In any case, the distribution differs significantly from a Bahcall-Wolf cusp, the 
predicted final state of a relaxed population. Given the age of the young giants, they must have formed 
rather close to the SMBH, either in-situ or transported in by in-falling clusters. Since the radial density 
distribution evolves only slowly on timescales of ~ few Gyrs this means that the young giants still contain 
information on their initial distribution. The 3D velocity of the young giants compared with the cool giants 
in the same magnitude bin is shown in Fig. 1121 The 3D velocity was computed, assuming a distance of 



8.3 kpc to the Galactic Center. The maximum 3D velocity w max = y 2G ^f m for stars to be bound is indicated 
in the figure. The projected distances r yields a lower limit for the physical 3D distance R. The kinematics 
of the young population are consistent with the old stars. Using a Kolmogorow-Smirnov test to asses the 
likelyhood that the young giants 3D velocities can be drawn from the cold population returns a probability 
of 70%. The fact that both giant populations share the same radial distribution is surprising. The young 
giants (Mzams ~ 3M ) are three times more massive than their old counterparts with a lifetime too short 
to experience relaxation and especially to change their angular momentum significantly. However, it renders 
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a radially varying IMF unlikely. The latter is a valid claim although star formation history and IMF are 
highly degenerate, since one effect would need to cancel exactly the other to emulate a radially constant 
distribution of young-vs-old giants. 



5. Calculation of the star formation history 



To derive the star formation history, we compared the observed distribution of giants in the H-R diagram 
to model populations. We restrict our analysis to the ttik < 15 magnitude range with > 50% completeness. 
This left us with 450 giants for the fit. We repeated the calculation with a completeness limit of 80%. 
The results were similar, yet with larger uncertainties. Thus we are confident that no bias is introduced 
due t o the completeness limit c hosen. The synthetic model populations were created with the IAC-STAR 
code (jAparicio &: Gallartl 120041 ) . The code allows a sele ction of differ e nt ste llar evolution and bolometric 
cor rection librarie s for t he computation. We chose the iBertelli et al.l (| 19941 ) stellar evolution library and 
the iLeieune et al.l (|1997l ) bolometric correction library since these are the only libraries with evolutionary 
tracks of stars with masses > 10 M Q . The code requires several input parameters such as metallicity, age 
and slope of the IMF. Since a ge and metallicity introduce some degeneracy, we adopted the mean metallicity 
found by ICunha et al.l (|2007j) . Thus, we set up model populations with a metallicity of Z = 1.5 Zq. For 
the calculation of the star formation history, we applied a method similar to iBlum et al.l (|2003l) . We defined 
four age bins (50-200 Myrs, 0.2-1 Gyr, 1-5 Gyrs and 5-12 Gyrs) with constant star formation rate within each 
age bin. The bins were chosen such that the evolutionary tracks are distinct enough to be resolved with 
the given data quality. The isochrones for stars with ages > few Gyrs are very similar. Therefore, the older 
age bins were chosen to be wider than the younger age bins. The younger age bins were selected by visual 
comparison of the model populations with the data. Fig. ITOl shows the chosen age bins. To study the impact 
of the IMF, we created models with various IMF slopes between —2.7 < a < —0.45. Among the models, we 
included a Chabrier/Kroupa IMF (a = -2.3), a top-heavy IMF (a = -0.85) jManess et al.ll2007 ) and a flat 
IMF with a slope of a — —0.45 as it was found recently for the young stellar disk (jBartko et al.ll2010l ). We 



note, that the star formation could have proceeded in episodic bursts. Yet, given the quality of the data, it 
is not possible to distinguish between a burst and a continuous formation within an age bin. Therefore the 
derived star formation rate is an average across the age bin. 



5.1. Fitting procedure 



We added Gaussian noise to the synthetic model populations, representing the statistical errors of 
-Mbd and T e fj. Then we removed stars from the synthetic populations according to the estimated complete- 
ness as a function of K-band magnitude. The completeness is not a function of temperature, since we used a 
minimum SNR criterion. We then binned the data and the model populations into an H-R diagram with a 
binsize of 1.5 a of the typical errors. The data were then fit as a linear combinati on of the H-R diag rams of the 
4 model populations (4 age bins) . For the fit we used the IDL routine TNMIN (|Markwardtll2009f ) . As a min- 



imization parameter we used the Poisson maximum likelyhood parameter \\ — 2 2? m i — n i + n i m ( n i/ m i)j 
where m, is the num ber of stars predicted by the model and rii is the number of observed stars in the ith bin 
of the H-R diagram ( Mighell 19991 Dolphin 2002 ). The contribution of each age bin determines the relative 
star formation rate of that bin. The best-fitting star formation rate for various IMFs is given in Table[3] 
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Fig. 11. — Ratio of young-to-old giants with magnitudes between < Mboi < —4 as function of radius. The 
average ratio is indicated (dashed line). The errors are the Poisson errors of the young giants in the given 
radial bin. The ratio is consistent with being flat, indicating that the young giants exhibit the same radial 
distribution as the old giants. 




Fig. 12. — The 3D velocity of the warm giants (red) with magnitudes of < Mboi < —4 and the cold 
comparison sample (gray) is shown. The maximum allowed velocity for stars to be bound is indicated; 
SMBH alone (solid line); SMBH + cluster mass (dashed line). 
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5.1.1. Uncertainty and quality of the fit 

To assess the fit uncertainty we used a method generally referred to as bootstrapping. We constructed 
1000 H-R diagrams by drawing random stars out of the data. Each star was allowed to be drawn any 
number of times. To include the systematic uncertainty of the data and the theoretical isochrones, we 
added a random temperature offset with a Gaussian a of 100 K to the whole data set. The constructed H-R 
diagrams were fit again with the model populations. The scatter of the derived relative star formation rates 
represents the la uncertainty of the relative star formation rate (see Table[3]). Note, the star formation 
rate uncertainty is mainly driven by the systematic temperature uncertainty. To judge the quality of each 
model fit, we produced 1000 Monte Carlo data sets, each containing 450 stars drawn randomly from the 
respective model (not the original data). Calculating the variance in each bin resulted in a variance diagram. 
We defined the fit quality as x 2 = Si ~" , where erf is the variance in bin i, mi is again the number 
of model stars and rij the observed number of stars in bin i. The values for % 2 , the number of degrees 
of freedom and the corresponding probability that the data can be drawn from the model population are 
represented in TableO The numbers of degrees of freedom arc equal to the number of populated bins minus 
the number of free parameters. Fig. [13] illustrat es the residual s for various models. Each panel shows the 



residual model - data weighted by the variance ([Dolphin! 120021 ). The color coding indicates bins where the 



model predicts more stars than observed (bright) respectively where the model underpredicts the number of 
stars (dark) . The first panel of Fig. [13] shows the best- fit model, while the second panel shows a continuous 
star formation (both assume a Chabrier IMF). The continuous model predicts too few old (i.e. cold) stars 
compared to the observations, while it predicts too many young (i.e. warm) stars. The data is well fit by a 
mostly old (> 5 Gyrs) population with an admixture of recently formed stars (< few lOOMyrs) as in case 
of the best-fit model. The discrepancy between a continuous formation scenario and a time varying star 
formation increases in case of IMFs favoring the formation of massive stars (right two panels). 



6. Results 

The data fit best a normal Kroupa/Chabrier IMF (49% acceptance probability, see Tab. [3]) but either 
flatter or steeper IMF can be accommodated within the 2 sigma uncertainties. The distribution in the H-R 
diagram does not strongly constrain the IMF, in contrast to the diffuse light and dynamical mass discussed 
in Sections 16.2.11 and 16.2.21 This is not surprising, since the old giant branch contains stars with zero age 
main sequence (ZAMS) masses between lM a < M < 1.2 M©. Thus the mass interval is too small to be 
significantly affected by changes of the IMF slope. The slope can affect the abundance of young giants, since 
they cover the mass range 2.5 Mq < M < 6M . However their abundance depends sensitively on the star 
formation rate of the last few hundred Myrs. Furthermore, the bright and massive end of the population 
suffers from low number statistics. This makes statements on the IMF slope uncertain. In general, the IMF 
and the star formation rate as a function of time are largely degenerate. Given the coarse sampling of the 
age bins, the data is not sufficient to constrain the IMF by the fitting procedure itself. The fit only allows to 
derive relative star formation rates for each assumed IMF. However, it is possible to constrain the IMF using 
the measured dynamical mass. Since each model predicts a certain mass composition of the population, it is 
possible to distinguish between different models. In any case, no acceptable fit was achieved for continuous 
star formation scenarios, irrespective of the assumed IMF. 
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Fig. 13. — Residual diagrams of two best-fit models (with a normal Kroupa IMF, a = —2.3 and a flat IMF, 
a = —0.45) compared with models of continuous star formation and the same IMFs. Each panel shows the 
residual model - data H-R diagram weighted by the variance in each H-R bin. Lighter areas indicate the 
model overpredicting the number of stars in the corresponding bin compared to the actual data. Darker 
regions indicate bins where the model predicts too few stars compared to the data. It is apparent that 
the data is well matched by models allowing for a time varying star formation rate. The continuous star 
formation models predict too few old stars on the cold giant branch, while they overpredict the number of 
young and massive giants on the warm branch. The discrepancy is worse in case of the flat IMF. The best-fit 
models however reproduce both branches and show only little discrepancies. 
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Table 3: The star formation rate for each age bin and various model IMFs. 



A/Tnrlpl 


Age oin [yjyi\ 


TMF olnnp a 


v 2 /d n f (Wnh "l b 


Ahsnlnte SFR h fl -4 M^nr -1 ] 


CSFR 


Model 1 


12-5 


-2.7 (-1.3) 


122/Q4 (3%1 


3 


0.6 




^ 1 








<i 




1 9 
± - u .z 






1 


^ 
u.o 




u.z - u.uo 






1 9 


Q 
O 


Model 9 (Chahrier) 


12-5 


-9 3 (-1 3) 


105/104 (4Q%) 


(J 


6 




5-1 








0.5 




1 - 0.2 






0.5 


0.5 




0.2 - 0.05 






8 


2 


Model 3 


12 - 5 


-1.5 


157/106 (0.1%) 


13 


3 




5 - 1 






2 


1 




1 - 0.2 






1 


0.7 




0.2 - 0.05 






11 


3 


Model 4 (Top-heavy) 


12 - 5 


-0.85 


143/100 (0.3%) 


126 


27 




5 - 1 






24 


11 




1 - 0.2 






6 


3 




0.2 - 0.05 






50 


12 


Model 5 (Flat) 


12 - 5 


-0.45 


131/107 (6%) 


614 


124 




5 - 1 






97 


43 




1 - 0.2 






22 


10 




0.2 - 0.05 






145 


34 


Continuous 


12 - 0.05 


-2.3 (-1.3) 


178/105 (< 0.01%) 


2 




Continuous 


12 - 0.05 


-1.5 


181/111 (< 0.01%) 


6 




Continuous 


12 - 0.05 


-0.85 


176/116 (< 0.01%) 


34 




Continuous 


12 - 0.05 


-0.45 


211/115 (< 0.01%) 


113 





"The upper and lower mass cutoff is: 0.5 < M < 120. The models 1 and 2 extend to 0.1 < M < 120 with a natter slope between 
0.1 < M < 0.5. 

^Outlier corrected \ 2 an d degrees of freedom. Outliers are bins off by more than 10 a. The d.o.f. corresponds to the number of 
populated H-R diagram bins (ignoring significant outliers) minus the number of fit parameters. Model 1-5 have four independent 
fit parameters (age bins), while the continuous models have only one free scaling parameter. In brackets, the acceptance 
probability is stated. 
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6.1. Star formation rate over cosmic time 

In the following we only consider the formation scenario with a Chabrier/Kroupa IMF. This is mo- 
tivated by the slight preference of the fitting, the intriguing total mass prediction of the model matching 
the dynamical measurements and the agreement with the observed diffuse background light (see next sec- 
tions). We find that the giant population of the Galactic Center is old. Fig. [14] shows the star formation 
rate as function of time (assuming a normal IMF). The formation rate more than 5 Gyrs ago was on average 
3 ± 0.6 x 10~ 4 M©/yr. According to our fit, roughly 80% of the total mass formed more than 5 Gyrs ago. 
The formation period was followed by a period of reduced star formation lasting another 4-5 Gyrs. The star 
formation rate reached a minimum about 1 Gyr ago. During the last few hundred Myrs, the formation rate 
increased again. The disk of young stars a re part of th e incre ased formation rate. Our results are in excel- 



lent agreement with the earlier findings of iBlum et al.l (|2003[) (compare Fig [T4l . The recent period of star 



formation accounts for about 10% of the total formed mass. Note, the recent star formation rate seems to 
be higher than several Gyrs ago. However the bin widths are very different. Each bin represents an average 
formation rate. The actual formation several Gyrs ago can have happened on short bursts with significantly 
higher star formation rates. The present day cluster contains roughly 50% of the total processed gas mass in 
living stars and 10% in remnants. The rest has been lost via stellar winds, explosions and potentially been 
swallowed by the SMBH. For simplicity, we set up a simple analytic model that approximates the measured 
SFR as function of time (solid line in Fig. [14]). 



SFR(t) = 6.8 x 10- 5 M Q /yr x e t/5 - 5 Gyrs + 4.3 x 10- 3 M Q /yr x e -*/o.06Gyrs_ ^ 

where t is the look-back time. The integration of the model yields the integrated mass as function of time. 
As Fig. [15] shows, about half of the cluster mass formed before a redshift of one. This suggests, that the 
nuclear cluster formed at a time, when the galaxies build up most of their stellar mas s. In that sense, th e 



nuclear cluster formed at the same time, or shortly after the bulge, about 10 Gyrs ago (|Zoccali et al.ll200 



6.1.1. Implications from the faint main-sequence population 

The two confirmed A-stars and the one candidate (Sec. l3.lTTj) are ideal tracers for star formation during 
the last few hundred Myrs. Early-type stars with K-band magnitudes 17 < rriK < 18 (assuming Ak = 2.8 
and Rq = 8.3 kpc) are main sequence B9/A0 dwarfs with average lifetimes of 500 Myrs. Late-type stars 
in the same magnitude range can only be giants that have already left the main sequence. These giants 
have ages between 1 to 12 Gyrs. Thus the number count ratio of the two populations provides a measure 
of the star formation efficiency during the last 500 Myrs relative to earlier star formation. This ratio is 
largely independent of systematic uncertainties since it avoids issues like incompleteness and spectroscopic 
dctectability. We included the candidate star, in the early-type population. In total we found three (two 
confirmed) early-type stars and 30 late-type stars in this magnitude bin. We used the IAC-STAR code and 
found that in case of continuous star formation over 12 Gyrs with a Kroupa IMF, the predicted early-to-late 
ratio is close to unity. The observed ratio therefore argues for a average star formation rate during the 
last 500 Myrs of only about 10 ± 6% of the average formation rate 1-12 Gyrs ago. To convert the relative 
formation rate into an absolute rate, we used the star formation rate inferred in the previous section and 
calculated an average rate of 2.3 ± 0.5 x 10~ 4 M© between 1-12 Gyrs look-back time. Using this value, 
we obtained an average star formation rate of 2.3 ± 2.1 x 10~ 5 M© during the last 500 Myrs (see Fig. lM]) . 
Within the errors this is consistent with the one derived by the H-R fitting. The A-star ratio provides an 
independent measurement of the relative star formation rates. It is another indicator for the early formation 
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Fig. 14. — Left: Star formation rate of the Galactic Center as function of time. The black circles represent 
the best-fit to the H-R diagram with a Chabrier/Kroupa IMF. The mass error is given by the la error of 
th e fit. The age error is simply the width of the age bin. For comparison, the star formation history derived 
by iBlum et al.1 (|2003l ) scal ed to a radius of 1 .2 pc is indicated (green triangles). The star formation rate 
derived by red supergiants (IGenzel et al.ll2010l) within 1 pc (red square, filled) and within 2.5 pc (red square) 
and the stellar disks (blue circle) (jBartko et al.ll2009l ) are indicated. The star symbol (purple) indicates the 
recent star formation rate inferred by the early-to-late ratio of faint GC stars fSec. l6.1.Tj) . The solid gray 
line shows a simple exponential model (see text) of the star formation as function of time. Right: The total 
mass formed in each age bin is shown. Although star formation occured at a high rate during the last few 
hundred Myrs, the total mass contribution is < 10%. The bulk of the stellar mass formed more than 5 Gyrs 
ago. 
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Fig. 15. — The integrated mass as function of time as obtained by Eq. 3 in Sec. 16. II is shown. The dashed 
curves represent the 1 a uncertainty. Note, today's cluster mass is only about 50% of the total formed mass. 
The remaining mass is lost due to stellar evolution. 



of the nuclear cluster. 



6.1.2. Comparison with previous work 



Maness et al.l (|2007l ) found that the giant population of the GC is on average warm and thus young. 



Consequently, they favored mod els with a normal IM F and an increasing star formation rate or a top-heavy 
IMF and continuous formation. iManess et al.l (|2007l ) used a CO index that is very sensitive to systematic 
effects. In particular, their CO index, t hough widely used, was recently discovered to vary with spectral 
resolution (|Marmol-Queralt6 et al.l 120081 ) . Since no suitable library was publicly available at the time of 
writing, they also used template libraries with resolutions between R ~ 3000 — 5000 to calibrate their 
R ~ 2000 data. The sum of the effects caused an underestimation of the CO equivalent width and an 
according overestimation of the stellar temper atures. This led to th e interpretation of a mostly young giant 
population. Applying the same methods as IManess et al.l (|2007l ) tog ether with a CO in dex that is less 
susceptible to systematic effects and using the newly available library of lRavner et al.l (|2009l ) with R ~ 2000, 
allowed us to revisit the star formation history of the Galactic Center. 



6.2. The mass composition of the nuclear cluster 

The present-day mass composition of the nuclear star cluster depends sensitively on the formation 
history and the IMF. The mass contribution of stellar remnants and stars as well as the total amount of 
consumed gas depends on the age and the IMF of a population. The consumed gas mass is larger than the 
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sum of remnant and stellar mass because a significant mass fraction is lost due to stellar evolution. Since 
massive stars lose a larger fraction of their initial mass, the gas consumption increases with flatter IMFs 
for the same total stellar mass. The same is true for the remnant mass. Massive stars have main-sequence 
lifetimes that are significantly shorter than the age of the galaxy. Therefore many generations of massive 
stars evolve through time and finally end as remnants. Thus, IMFs that favor massive stars produce more 
and heavier remnants. As a consequence, the mass contribution of stars still burning hydrogen drops with 
flatter IMF slopes. 



6.2.1. Constraints from the dynamical mass 



Various models have been fit to the data (Sec. l5.ip . Each model yields the star formation rate as function 
of time under the assumption of a certain IMF. To derive the absolute mass contribution, it is necessary to 
scale the models according to the actua l star counts. To get a representative number for the nuclear cluster 
we used the stellar surface density from lSchodel et al.l (|2007l ). As discussed in Sec. Q] the radial extent of the 
NSC is ~ 5pc. Our data, however, probed only the inner 30" (1.2 pc). Therefore we restricted ourselves to 
that radius. By assuming 30" as a sharp edge, the surface density from Schodel et al. (2007, their fig. 12) 
yields ~ 18200 stars with hik < 17.75 within a projected radius of 30" from Sgr A*. Between 45% to 55% 
of the proj ected stars are also co ntained within a 3D distance of 30", depending on the radial density profile 
derived by ISchodel et al] (|2007l) . Thus, we assumed a total of 9100 stars with R 3 d < 1.2 pc to scale the 
remnant-, stellar- and total gas mass for each model. Table[3] shows the derived star formation rate for each 
model IMF. Fig. [16] shows the mass composition for the models derived in Sec. 15.11 The model predictions 
have to be compared with the dynamical mass estimates of the nuclear cluster (see Fig. [T6|) . The total 
mass enclosed within 30" is 5.7 ± 0.9 x 10 6 M . The SMBH contri butes 4.3 ± 0.5 x 10 6 M^ dGhez et al 
20081: iGillessen et al.ll2009h. while the remaining 1.4 ± 0.7 x 10 6 M (jGenzel et al.lll996t iTrippe et alJl200 



Schodel et al 



2009UGenzel et alJboiol) are contained in stars and stellar remnants. Numerical simulations of 



Freitag et al. ( 20061 ) have shown that during 10 Gyrs about 4 x 10 5 M stellar and 1 x 10 5 M remnant mass 
is removed from the cluster by the SMBH through tidal disruption and inspiral. The cluster composition 
is scaled to match the observed number of giants. Thus, the inferred stellar mass intrinsically accounts for 
stars lost to the SMBH because giants and giant progenitors are equally likely to be disrupted as unresolved 
main-sequence stars constituting the bulk of stellar mass. The situation for remnants, however, is different. 
The main remnant contribution comes from stellar BH dominating the mass density in the inner 0.1 pc. Due 
to their mass (~ 1OM ) they are significantly affected by mass segregation and more likely to inspiral into the 
SMBH. Consequently, the inspiral of stellar BHs has to be taken into account in the mass budget. We simply 
added 10 5 M with an u ncertainty factor two to the observed dynamical mass to constrain the IMF. We took 
the value determined bv lFreitag et al. (2006) as being constant with IMF. This might be the weak point in 
the line of argument, since iFreitag et al.l (|2006f ) assumed a Chabrier/Kroupa IMF in their simulations. As 
Fig. [16] shows, the cluster composition depends sensitively on the assumed IMF. For flatter IMFs, the number 
of stellar BH increases, however, the mean mass increases accordingly making mass segregation less efficient. 
Thus the actual mass transfer between the cluster and the SMBH as function of IMF can only be addressed 
by further simulations. Keeping the limitations in mind, the dynamical mass plus the removed mass therefore 
make IMFs with a slope flatter than a = —1.1 unlikely. Even if the transferred mass is significantly higher, 
models flatter than a = —0.8 violate the total enclosed mass (SMBH+cluster). Mechanisms that expel 
remnants or stars are very inefficient. One mechanism is the evapor ation of stars fro m the nuclear cluster. 
However, for the GC the timescale is greater than a Hubble time ( Alexander 20051) . Mass segregation is 
another mechanism for removing stars. This mechanism moves low-mass constituents outward, and high- 
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mass c onstituents inward, without changing the radial density profile significantly. iMiralda-Escude &: Gould 
(|2000h suggest that inward migration of a considerable number of bulge stellar blackholes can remove low 
mass stars from the nuclear cluster. According to their model, the low mass stars are supposed to form a core 
with a radius of 1-2 pc. Observations however find a core radius of about 0.3 pc (jSchodel et al.ll2007l) . If the 
observed core is related to mass segregation, then the effect is smaller than predicted. Our constrains on the 
IMF are however not weakened by potential mass segregation. Our models use the number of observed stars 
as a scaling for the total processed mass (and remnants produced) . Therefore if low mass stars are removed 
from the inner region, our scaling underestimates the remnant mass. The inferred stellar mass however is 
not changed, because giants and unresolved main sequence stars have very similar masses (~ 1M Q ). Thus 
the total predicted mass will grow if we take into account that low mass stars are removed due to mass 
segregation. This makes our statement on the IMF constrains even more significant. In general, it is very 
hard to remove mass in form of stars and remnants from a nuclear cluster (a SMBH binary, however is 
able to remove stars efficiently). The processed gas mass, unlike stars and remnants, can be expelled from 
the nuclear cluster due to supernovae or AGN activity in the distant past of the SMBH. Stellar winds are 
also able to remove gas from the cluster. While the low mass main-sequence stars and giants have wind 
speeds significantly lower (u w ind ~ few lOkm/s) than the escape velocity of the NSC the most massive 



stars (fwind ~ 1000 km/s) can efficiently remove gas from the cluster ([McLaughlin et al 



2006; 



Martins et al 



20071 ). Therefore, the processed gas mass cannot constrain the IMF although flat IMFs require a factor 



hundred greater gas masses than Kroupa like IMFs. Summarizing the arguments we can say that IMFs 
flatter than a > —1.1 violate the observed cluster mass. IMFs flatter than a > —0.8 exceed even the total 
mass (SMBH+cluster). 



6.2.2. Constraints from the diffuse H-band background 

As discussed in Sec. l2.1.fl depending on the IMF and star formation history the diffuse background light 
from faint unresolved stars may probe the bulk of the stellar population. We find that the diffuse background 
contributes H^g/H — 27±9% to the total H-band flux. Furthermore, we derived a mass-to-diffuse-light ratio 
of Mdyu/i/diff = 2.6 ± 1.5 Mq/ Lh q. In the following we used the inverse ratio -ffdiff/Miyn f° r convenience. 
We compared both quantities to our best-fit models fSec. l5.1[) . To show how H^g/H and H^g/M depend on 



the I MF slope and overall star-formation history, we used the population synthesis code STARS ^Sternberg 



1998) to compute these ratios for a range of IMFs and star- formation time-scales to for simple histories 
SFR(t) oc e~*/*°. We considered a Kroupa IMF, and also three power-law IMFs, dN/dm cx m~ a , with a 
equal to -1.5, -1.35, and -0.85. All of the IMFs range from 0.01 to 120 M Q . We used Geneva evolutionary 
tracks for solar metallicity stars, combined with empirical colors and bolometric corrections for dwarfs, giants, 
and supergiant stars. For low-mass stars (< 0.8 M Q we computed the H-band luminosities along the lower 
main-sequence using the calibrations of Henry & McCarthy (1993). We considered exponentially decaying 
(t = 3 Gyr), continuous (t — oo) and exponentially increasing (t = —3 Gyr) star-formation rates, for an 
assumed cluster age t = 13 Gyr, and fixed dynamical mass Md yn — 1.5 x 10 6 M Q . The results are displayed 
in Figure [T71 in the H^g/M versus H^g/H plane. The gray areas indicate the measurement including the 
1 er uncertainty. Figure [T71 shows that H^g/M and H^g/H both decrease for flatter IMFs. This behavior is 
due to the relative increase in remnant mass from massive stars, and the reduction in the relative fraction of 
diffuse light from low mass stars. For a given IMF, H^g/M decreases but H^g/H increases, as the history 
is altered from increasing, to steady, to declining star-formation, because the total luminosity for fixed mass 
is reduced for this sequence of histories, while the fraction of light produced by the accumulating long-lived 
low- mass stars is increased. In FigurefTTl we also indicate the positions of the "best-fit" models discussed in 
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Fig. 16. — Mass composition of the inner 1.2 pc after 12Gyrs of star formation. The composition and 
total mass depends sensitively on the assumed IMF. The symbols show the five fitted models (Tab. [3]) with 
various IMF slopes —2.7 < a < —0.45. Green triangles represent the stellar mass, blue squares represent 
the remnant mass, red open circles show the sum of remnant and stellar mass and black filled circles show 
the total processed gas mass. The dark gray horizontal line indicates the dynamical mass within 30 "and its 
uncertainty (light gray area). 
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Sec. 15.11 and listed in Tabled FigurefTTl shows that the observed H^g/M and Hdm/H ratios require steep 
(canonical) IMFs and continuous or mildly decaying star-formation. In particular, fiat IMFs with a > —1.5 
are inconsistent with the observations 



Our measurements of the diffuse H-band background can be compared to the findings of lLockmann et al 



(|2010l ). who used the diffuse K-band background together with the dynamical mass to derive M/K^m — 
1-Cto 7 -^©/-k© ,Ki which was then compared with model predictions. Their mass-to-diffusc-light ratio is 
somewhat lower than our value (M/K^is = 1.9 ± 1.2 M^/L^.x ), yet consiste nt within the errors. The 
di fference originat e s from the assumed extinction value. Lockmann et al. ( 201ot ) used Ak — 3.3 published 



bv lBuchholz et al. ( 2009). Howev e r, this value is 0.5 magnitudes (factor 1.6) higher than the most recent one 



published by Schodel et al. ( 2010l ): Fritz et al. ( 2011 ). Correcting for this difference removes the discrepancy 
between their finding and ours. Their analysis favored an IMF steeper than a = —1.3 together with a 
continuous star formation rate or an increasing star formation rate. Our analysis confirms their finding of a 
steep IMF. Taking i nto account the r e vised extinction our finding of an early star formation is also consistent 
with the findings of lLockmann et al] (|2010l ). However, we additionally used the ratio of diffuse-light to total 
light to constrain the star formation history even further. 



7. Discussion 

All statements concerning the star formation history of the NSC have been made under the assumption 
that the traced population is representative for the whole cluster. The radial extent of the cluster is ~ 5pc, 
i.e. larger than the radii probed with our SINFONI observations. However, the inner 1.2 pc contain already 
6 x 10 6 Mq of the total 30 x 10 6 M Q . Thus we cover about 20% of the total cluster mass and gives us 
confidence that our results are somewhat representative for the whole cluster. 



7.1. Initial Mass Function in the GC 



We find that the old giant population must have formed with a IMF steeper than a < —1.3, probably 
with a normal Chabrier/Kroupa IMF. In any case, the IMF mu s t have been significa ntly different from the 
IMF observed in the young stellar disks. Paumard et al. ( 2006 ); Bartko et al. ( 2010l ) found that the young 
stellar disks in the GC have formed p referentially massive stars. They favored an almost fiat IMF with 
dN/dm ~ to" - 45 (jBartko et al.ll2010| ). This indicates that the environmental conditions at the time of 



formation were different in the past. The nuclear cluster probably formed at a time, when the SMBH itself 
was younger and less massive. With the radius of influence smaller than today, the gravitational potential 
resembled more a normal cluster. Consequently the IMF was closer to the universal one. The nuclear cluster 
might also have been contaminated by in-spiraling clusters that had formed outside of the Galactic Center. 
Depending on their mass, clust ers within several 10-100 pc can spiral into the NSC during a Hubble time 
(jAgarwal fc Milosavlievial2011l ). Both effects lead to an IMF significantly steeper than the one observed in 
the young stellar disks today. 
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Fig. 17. — The observed ratios of diffuse background light to the total H-band luminosity i/diff /H (horizontal 
bar shows 1 a confidence interval) and the inverse mass-to-diffuse-light ratio -f/diff /Mdyn (vertical bar shows 
f a confidence interval) constrain the IMF and the star formation history. The colors indicate the assumed 
IMF. The best-fit models (Tab[3]) are represented by open circles. For comparison, simple models with 
exponentially declining, constant and exponentially rising (connected by solid, dashed, dotted lines) star 
formation are shown. The star formation rate in this case is parametrized as SFR(i) oc e _t /*°, where to = inf 
corresponds to continuous star formation, to = 3 Gyr corresponds to a burst with an exponential decline 
timescale of 3 Gyrs and to = —3 Gyr corresponds to a raising star formation rate. The assumed cluster 
age is t = 13 Gyr. Models with an increasing star formation rate fail to produce the observed amount 
of background light, while models with decaying star formation underpredict the fldiff /-Wdyn ratio. The 
mismatch gets worse for flatter IMFs. Only the best-fit or the continuous model with a normal Kroupa IMF 
match both constraints. 
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7.2. Star formation in the vicinity of the GC 



The Galactic Center region shows star formation at all ages. 



clusters at distances of a few ten pc were formed 2 and 4 Myrs ago (IFiger et al 



The massive Arches an d Quintuplet 

The stellar 



disk in the central cluster was formed about 6 Myrs ago (jPaumard et al 



200 



6; 



2002 



1999) 



Bartko et al 



20091) and at 



rough ly 50 pc distance the SgrB2 cloud harbors several massive star cluster in the making. iFiger et al 
(|2004l ) used HST photometry of several fields withi n 100 pc of the cente r to derive the star formati on history. 
They find a continuous star formation in their FoV. lGenzel et al] (|2003f ) and lBuchholz et al.l ([20091 ) used AO- 
assisted photometry of the nuclear cluster and found that the KLF is well matche d by a bulge-like (8-10 Gyr) 
population with an admixture of young main-sequence stars. iBlum et al. I (|2003h used spectroscopic data of 
the most luminous giants within 2.5 pc from SgrA*. Our work confirms their favored formation scenario. 
The star formation happened predominantly at old times. An intermediate period showed a reduced star 
formation, while during the last few 100 Myrs the star formation rate was increasing. This is also consistent 
with findings based on the stellar disks and the number of red supergiants. Roughly 75-90% of the mass 
contained in the central cluster formed 5-12 Gyrs ag o. The formation of the cluster might coincide with 
the formation of the bulge around 10 ± 2.5 Gyrs ago ( Zoccali et al. 20081 ). The a-element ratios require the 
bulge to have formed on a short timescale (~ lGyr). Whether this is also true for the Galactic Center 



cannot be answered with the given uncertainty of the data. The fact that iBuchholz et al.l ([20091 ) do not find 
a young giant population is not surprising, since those stars make up only ~ 10% of the number counts. 
The KLF fitting is not sensitive enough and since the young giants are rather warm and bright, they might 
be attributed t o the admixture of young main-sequence stars. The same is true if the results are compared 
\2Q0A 



to 



Figer et al 



They, however, find a continuous star formation rate within 100 pc. This does not 
contradict our findings is consistent wit h the idea of recurrent episodic star formation in the greater Galactic 
Center region ( Serabvn fc Morris] 1996h . The star formation events themselves happen on spatial scales of 
< 1 pc (e.g. Central-, Arches & Quintuplet Cluster). If one assumes that the individual star formation 
events happen incoherently , then one would naturally expect a continuous star formation history on scales 
>> 1 pc. Thus lFiger et al. (2004) find a star formation history that is smoothed across many episodic single 



events. 



8. Conclusion 



We used several methods to constrain the star formation history of the Milky Way Nuclear Cluster. This 
is the only nuclear cluster in which individual stars can be resolved and reliable age estimates can be made. 
For this purpose, we used 450 K-band spectra of late-type giants to derive individual stellar temperatures. 
By using a CO index that is insensiti ve to systematic eff ects such as reddening or instrument resolution, 
together with the new stellar library of lRavner et al] ([20091 ) we improved the temperature calibration for the 
red giants in the Galactic Center. Together with K-band photometry we were able to construct a detailed 
H-R diagram of the giant population. The comparison of the observed H-R diagram with model populations, 
allowed us to infer the star formation history of the Galactic Center. Our results are as follows: 



1. The bulk of the Nuclear Cluster is old. Roughly 80% of the stellar mass formed more than 5 Gyrs ago. 
It might have formed at the same time as the Galactic Bulge at a redshift of 1-2. 

2. After the bulk of the cluster had formed, a period of reduced star formation followed between 1-5 Gyrs 
ago. The star for mation reach e d a m inimum ~ 1 Gyr ago. Our inferred star formation history confirms 
earlier findings of Blum et al. ( 20031 ) . 
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Only during the last 200-300 Myrs star formation has set in at a significant level. A population of 
intermediate age giants are tracers of that period. Making up only 10% of the number counts, the 
intermediate age population is hardly traceable with the KLF. However, they can be clearly identified 
in a H-R diagram due to their high temperatures. The spatial distribution and kinematics of those 
giants resemble the one of the old giant population. 

We report the first detection of main-sequence B9/A0 stars with magnitudes 17 < rriK < 18 in the 
GC. They are the faintest early-type stars having been found in the Galactic Center so far. With these 
stars we probe the mass regime between 2.2 — 2.8 Mq and main-sequence lifetimes of 360-730 Myrs. 

The ratio of late-type to early-type stars in the magnitude bin 17 < juk < 18 yields an independent 
estimate of the star formation rate during the last 500 Myrs. We find that the average rate during the 
last 500 Myrs must have been a factor 10 lower than the average rate between 1-12 Gyrs ago. This 
finding supports the claim of an old cluster population. 

We find that the bulk of the stellar mass must have formed with an IMF steeper than dN/dm ~ 
m a ; a < —1.5. Otherwise, the required remnant and stellar mass violates the observed dynamical 
mass and diffuse background. Thus, the bulk of the old stars formed with an IMF significantly steeper 
than the one observed in the young stellar disk. We suggest that this apparent discrepancy can 
be naturally explained if the stars formed at a time when the SMBH itself was much younger and 
less massive. Consequently the sphere of influence was significantly smaller. Without the extreme 
environment of a SMBH, the stars fo rmed with an IMF close to the universal one. The SMBH was 
fed by gas including stellar mass loss ( Freitag et al. 2006) and dominated the inner few pc over time. 



With the growing sphere of influence, the IMF became flatter and reached the value observed in the 
young disk 

The deep census of the Galactic Center using integral field spectroscopy yields a late-type KLF between 
12 < rriK < 1 8 with a slope of d log-ZV/dw ^ = 0.33 ± 0.03. This confirms the previous findings 
of"" 



Genzel et al 



|Alexanderlll999h 



OK 

(|2003[ ); iBuchholz et all (|2009l) . The slope is consistent with an old Bulge like KLF 



Taking into account systematic effects leads to an improved age estimate of the Nuclear Cluster. To improve 
the age estimate further, however, requires a technical leap. The main difficulty is that > 5 Gyrs old 
isochrones are spaced by only a few 10 K. The necessary temperature accuracy can hardly be achieved 
by spectral analysis. The most promising way is to detect a turnoff of the main-sequence population. This 
requires the spectroscopic identification of stars with magnitudes of rriK ~ 19. For example, the confirmation 
of late-type giants only 1 mag fainter than the current limit can unambiguously confirm the presence of stars 
with ages of > 10 Gyrs. Yet, the required telescope time, which is needed for the spectroscopic identification 
of those stars, is tremendous. Furthermore, these faint stars are significantly affected by stellar crowding. 
Thus further progress will require the resolving power of a new generation of telescopes like the E-ELT or 
TMT. The same is true to derive a model-free measurement of the initial mass function in the vicinity of 
the SMBH. 



This work has made use of the IAC-STAR Synthetic CMD computation code. IAC-STAR is supported 
and maintained by the computer division of the Instituto de Astrofisica de Canarias. AS thanks the DFG 
for support via German-Israeli Project Cooperation grant STE1869/1-1.GE625/15-1. We thank the referee 
for his valuable comments. 
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